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ABSTRACT 

^ ; A scenario is explored in which Lya absorbers at low redshift arise 

from lines of sight through extended galaxy disks, including those of 



> 



dwarf and low surface brightness galaxies. A population of galaxies is 

simulated based upon observed distributions of galaxy properties, and 

^ ; the gas disks are modeled using pressure and gravity confinement. Some 

O ! parameter values are ruled out by comparing simulation results with the 

Q\ '. observed galaxy luminosity function, and constraints may be made on 

^ ■ the absorbing cross sections of galaxies. Simulation results indicate that 

. . . . . 

^ ; it is difficult to match absorbers with particular galaxies observationally 

d from luminous galaxies. Low impact parameter absorption is dominated 

by low luminosity dwarfs. A large fraction of absorption lines is found 
^ ' to originate from low surface brightness galaxies, so that the absorbing 

5^ \ galaxy is likely to be misidentified. Low redshift Lya absorber counts 

can easily be explained by moderately extended galaxy disks when low 
surface brightness galaxies are included, and it is easily possible to find 
a scenario which is consistent with observed the galaxy luminosity func- 
tion, with low redshift Lyman limit absorber counts, and with standard 
nucleosynthesis predictions of the baryon density, il^. 



Subject headings: galaxies: fundamental parameters, luminosity function- 
intergalactic medium — quasars: absorption lines 



1 



LiNDER 



Comparing Galaxies and Lva Absorbers 2 



1. Introduction 

The numerous absorption lines seen short- 
ward of hydrogen Lya in quasar spectra are 
a powerful probe of gas in galaxies. Recent 
observations of the Lya forest at low red- 
shift (Bahcall et al. 1996) have allowed for 
studies which compare absorbers with nearby 
galaxies, but the nature of the lower column 
density Lya absorbers remains controversial. 
Here I explore the possible connection of Ly- 
man limit systems {Nhi = 10^'^''^ to 10^°-^ 
cm~^) and Lya forest absorbers {Nhi < 10^^'^ 
cm~^) with galaxies at low redshifts by using 
simulations to relate properties of the galaxy 
population to absorbers. A complete under- 
standing of galaxy formation and evolution 
will require knowledge of the properties of 
the low redshift galaxy population, includ- 
ing galaxies which are more difficult to de- 
tect. Absorbing galaxies should include dwarf 
and low surface brightness (LSB) galaxies. 
These objects are generally found to be rich in 
gas, so they must make at least some contri- 
bution to Lya absorption. Numerous stud- 
ies have found extended gaseous envelopes 
around dwarf galaxies, so that their absorb- 
ing cross sections might not be small as their 
name would imply. 

Burbidge et al. (1977) first noted that 
known galaxies must have huge cross sec- 
tions in order to explain the Lya absorbers. 
However, galaxy disks may be quite extended 
when ionized gas is considered, and the ab- 
sorption capabilities of the more recently dis- 
covered dwarf and LSB galaxies have not been 
explored. Recent studies have had some suc- 
cess in matching absorption lines with partic- 
ular galaxies, although conclusions vary de- 
pending upon the column density threshold. 
Studies concentrating on high column density, 
low redshift absorbers generally find clear 



associations with galaxies (Lanzetta et al. 
1995a; Stocke et al. 1995; Steidel et al. 1997), 
while lower column density absorbers tend to 
be located far from any observed galaxy. Of- 
ten one or more galaxies will be seen at several 
hundreds of kpc from the line of sight, so that 
it is unclear whether any of the galaxies are re- 
sponsible for the absorption, although the ab- 
sorbers are generally found to follow the large 
scale galaxy distribution (Le Brun, Bergeron, 
& Boisse 1996; Morris et al. 1993; Bowen, 
Blades, & Pcttini 1996; Ranch, Weymann, & 
Morris 1996; Mo & Morris 1994; Salzer 1992; 
van Gorkom et al. 1996). A few absorbers 
have been found in cosmic voids (Stocke et 
al. 1995; ShuU, Stocke, & Penton 1996), while 
limited attempts have been made to detect 
LSB galaxies near quasar lines of sight (van 
Gorkom et al. 1996). Several groups have 
tested for a physical connection between ab- 
sorbers and galaxies by looking for an anti- 
correlation between the equivalent widths of 
the absorption lines and the impact param- 
eters between the galaxies and lines of sight 
(Lanzetta et al. 1995a; Bowen et al. 1996; 
Le Brun et al. 1996). However, it is not pos- 
sible to know that any particular absorption 
line has been matched to the correct galaxy, 
and some galaxies which are unidentified may 
be causing absorption. Thus it is also useful 
to estimate how many absorbers could arise 
from lines of sight through galaxies. 

Many nearby galaxies have been discovered 
only recently due to a severe selection effect 
against LSB galaxies (Bothun, Impey & Mc- 
Gaugh 1997). Disney (1976) made the first 
quantitative analysis of this selection effect 
after realizing that galaxies obeying the Free- 
man (1970) law were within one magnitude of 
the typical night sky brightness. The primary 
goal of recent surveys for LSB galaxies has 
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been simply to obtain a representative sample 
of the nearby galaxy population (Schombcrt 
et al. 1992; Impey et al. 1996; Sprayberry 
et al. 1996; Dalcanton et al. 1997a; O'Neil, 
Bothun, & Cornell 1997). Most LSB galax- 
ies exhibit a wide range of properties such as 
sizes and morphologies which are similar to 
those of high surface brightness (HSB) galax- 
ies. The exceptions are Malin objects, which 
make up about 10% of the Schombert et al. 
(1992) catalog. These objects are highly ex- 
tended and extremely low in surface bright- 
ness, and they are often among the most lu- 
minous galaxies known. While most surveys 
for LSB galaxies have been biased toward blue 
objects, varying colors have been observed by 
O'Neil etal. (1997). McCaugh (1996a) found 
that the distribution of central galaxy sur- 
face brightnesses is approximately flat to at 
least jiQ — B mag arcsec"^. Evidence has 
been seen more recently for a falloff faintward 
of 24 B mag arcsec'^ (O'Neil ct al. 1997; 
Zwaan, Briggs, & Sprayberry 1997; Briggs 
1997). O'Neil et al. (1997) attribute this ef- 
fect to an inability of extremely low surface 
brightness galaxies to form in a cluster envi- 
ronment. 

The gaseous extent of galaxies is highly un- 
certain. Several studies involving 21 cm map- 
ping of galaxy disks have found 'sharp edges,' 
where the neutral hydrogen column density 
falls off quickly from a few times 10^^ cm~^ 
to an undetectable level. Such 'edges' have 
been explained successfully by models where 
the ionization level increases rapidly from the 
optically thick to optically thin regime (Mal- 
oney 1993; CorbeUi & Salpeter 1993; Dove 
& ShuU 1994), thus making a sudden cutoff 
in the total column density of the gas un- 
necessary. A detection of ionized gas beyond 
the HI disk in NGC 253 has been made re- 



cently (Bland- Hawthorn, Freeman, & Quinn 
1997). Ionized outer disks could be quite ex- 
tended if they arc not severely affected by 
galaxy interactions. These extended struc- 
tures may also be consistent with results of 
cosmological simulations, where galaxies form 
within extended sheets and filaments (Wein- 
berg, Katz, & Hcrnquist 1997; Norman ct al. 
1997). Many absorbers are found to arise 
in the sheet structures which have already 
formed at high redshifts. Double hne of sight 
observations (Dinshaw et al. 1994, 1995; Fang 
et al. 1996) have shown that absorbers are 
hundreds of kpc in size. Furthermore, these 
observations were found to be consistent with 
flattened absorbers having organized motion, 
such as rotating disks (Charlton, Churchill, & 
hinder 1995). Disks in reality are likely to be 
warped and broken, but when considering the 
average behavior of a large sample of galaxies, 
some systematic falloff in total column den- 
sity with radius can be expected. Only fleld 
galaxies are considered here, without the ef- 
fects of clustering and interactions. 

The number density of galaxies at low red- 
shift is also uncertain due to a possible abun- 
dance of LSB galaxies. An upper limit can 
be found by supposing that all absorbers are 
caused by galaxy disks. The average num- 
ber of absorbers per unit redshift along a line 
of sight was found by Bahcall et al. (1996) 
to be {dN/dz)o — 24.3 ± 6.6, complete to an 
equivalent width in Lya of 0.24 A. The neu- 
tral hydrogen column density distribution is 
observed to be approximately a power law at 
high redshifts (but see Petitjean et al. (1993) 
and Charlton, Salpeter, & Linder (1994)) 
with a slope of -1.5. While an accurate distri- 
bution of neutral hydrogen column densities 
at low redshift has not been published yet, it 
is assumed to look similar to those at higher 



LiNDER 



Comparing Galaxies and Lva Absorbers 4 



redshifts. A significant change in slope is gen- 
erally not observed between redshifts 4 and 
2.7 for the column density range of interest 
here (Lu et al. 1996; Hu et al. 1995; Kirk- 
man & Tytler 1997; Kim et al. 1997). 

Section 2 of this paper describes the method, 
including defining and simulating galaxy pop- 
ulations and modeling the galaxy disks as ab- 
sorbers. Results are discussed and compared 
with local luminosity functions, the predicted 
baryon density, and absorber/galaxy observa- 
tions in section 3. Future papers will dis- 
cuss the effects of clustering and the possi- 
bility of distinguishing between this scenario 
and and others, such as those with significant 
nongalactic absorbers. The value of Hq — 100 
km s~^ Mpc~^ is assumed. 

2. Method 

Samples of galaxies were simulated from a 
galaxy population as defined in section 2.1, 
where the method for choosing galaxy lumi- 
nosities and surface brightnesses is explained. 
In each simulation, all of the galaxies were 
given random coordinates within a cube of 
space with a side length of 50 Mpc. Each 
galaxy disk was modeled using pressure and 
gravity confinement, as described in section 
2.2. Lines of sight, going straight through the 
cube, were then chosen randomly. When a 
line of sight intersected a galaxy disk, the neu- 
tral column density was found by integrating 
the neutral gas density along the line of sight. 

2.1. Defining a Galctxy Population 

A galaxy population is defined here us- 
ing simple parameter distributions which are 
consistent with observations. Galaxy lumi- 
nosities are determined by the choice of scale 
lengths and central surface brightnesses. Re- 



lations are then assumed between these ob- 
served parameters and the properties of the 
gas disks. 

The surface brightness distribution is as- 
sumed to be fiat with an exponential falloff 
at the high surface brightness (HSB) end, as 
in McGaugh (1996a). The LSB end is known 
to extend to at least 24 B mag arcsec"^, 
while a few objects have been observed with 
Ho >27 B mag arcsec"^. Surface brightness 
profiles for the individual galaxy disks are as- 
sumed to be exponential, where the distribu- 
tion of observed scale lengths is found to be 
a power law, (pihs) oc k^, where n ~ — 2 
(van der Kruit 1987; Hudson & Lynden-Bell 
1991; de Jong 1995). No correlation is as- 
sumed between scale length and central sur- 
face brightness, except that large, HSB ob- 
jects are known not to exist (Sprayberry et 
al. 1995). Bothun, Impey & Malin (1991) 
have shown that the above distributions re- 
quire a slope of —1.5 in the faint end of the 
galaxy luminosity function. 

The bright end of the luminosity function 
is generally well fit by a Shechter function 
(Marzke, Huchra, & Geller 1994). Undiscov- 
ered Mahn objects cause some incompleteness 
(McGaugh 1996a), although they are thought 
to be rare (Dalcanton, Spergel, & Summers 
1997b: Hoffman, Silk, & Wyse 1992). Addi- 
tional faint galaxies have been added in evolu- 
tionary models which avoid strong evolution 
in the faint blue galaxy population at higher 
redshifts (Gronwall & Koo 1995; McGaugh 
1994; Driver et al. 1995). Sprayberry et al. 
(1997) found that a combined Shechter bright 
end and power law faint end {Mb > —16) 
provide the best fit for the observed luminos- 
ity function, including LSB galaxies, where 
the faint end slope was found to be quite 
steep (~ —2). Observed luminosity functions 
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(Driver & Phillips 1996; Zucca et al. 1996; 

Marzke et al. 1994) are generally found to be 
consistent with such hybrid shapes, includ- 
ing a Shechter function with a = —1 and a 
power law with slope between —1.4 and —1.8 
for Mb > —16, as discussed in Sprayberry et 
al. (1997). Here the disk scale length dis- 
tribution, for a given fiQ, obeys a power law 
if Mb > Mbsw For more luminous galaxies, 
a scale length is chosen, for a given central 
surface brightness, so that the Marzke et al. 

(1994) Shechter function is obeyed. 

LSB galaxies obey the Tully-Fisher rela- 
tion (Zwaan et al. 1995), and simulations by 
Navarro (1997) have hinted at a physical ex- 
planation for this relation. Here the Tully- 
Fisher formula from Jacoby et al. (1992) 
is used to determine a rotation velocity for 
each galaxy, which is needed to model the 
dark matter halo. Evidence is given by de 
Blok, McGaugh, & van der Hulst (1996) that 
LSB galaxies are low in central column den- 
sity compared to HSB objects. A central 
column density value for HSB galaxies has 
been estimated by Bowen, Blades, & Pettini 

(1995) . A relationship is assumed here which 
is consistent with the Bowen et al. (1995) 
estimate and which requires a factor of five 
change in surface brightness for a factor of two 
in central column density (McGaugh 1996b) 
so that the central column density becomes 

2.2. Modeling GalcLxy Disk Absorbers 

The total (neutral plus ionized) column 
density is assumed to fall off exponentially in 
each inner disk and as a power law farther 
out: 



where r is the distance in the plane of the disk, 
is the value of r at which the switch from 
exponential to power law occurs, and /i2i is 
an exponential scale length of the inner neu- 
tral hydrogen disk which could be determined 
using 21 cm mapping. Evidence for a power 
law falloff in outer galaxy disks has been seen 
by Hoffman et al. (1993), and such behavior 
is required for galaxies to produce a power 
law type neutral column density distribution. 
A neutral column density distribution which 
goes as N^^'^ corresponds to p = 4. The value 
of Tgyj is quite uncertain, and here it is as- 
sumed to be at one or two times the radius 
at which the gas becomes highly ionized (r^). 
The ratio of the neutral hydrogen scale length 
to B scale length (/?.2i/^b) is also uncertain. 
An initial estimate of 3.4 is used here based 
upon the conjectures of Salpeter (1995), and 
a smaller value (1.7) is also tested. 

The disks are modeled using pressure and 
gravity confinement (Charlton, Salpeter, & 
Hogan 1993; Charlton et al. 1994). The ver- 
tical structure is calculated in each disk, so 
that the internal pressure and density are re- 
lated by P{z,r) = 2ntot{z,r)kT, where the 
assumed temperature is T = 20000 K, z is 
the height above the disk plane, and Utot is 
the density of neutral plus ionized hydrogen. 
The inner disks are confined by a spherical 
dark halo with gravitational acceleration g = 
V^gfZ/r'^ where Vrot is the rotation velocity. 
The boundary condition P(z — w/2) = P^xt 
defines the absorber width w. The outer 
disks are confined by an external pressure P^xt 
which may be due to infalling material or to 
a local or diffuse intergalactic medium. The 
high redshift value of Pext/k — 10 cm~^ K 



Ntot{r) 



iVoexp(-r//i2i) for r < r^^ 

NQeyiY>{-rsw/h2i){r/rswY^ for r > Vgw 

(1) 
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(Charlton et al. 1994) is used here, al- 
though more recently evidence has emerged 
supporting a similar value at low redshift 
(Wang & Ye 1997). Hydrostatic equilibrium 
is assumed, so that 



dP 

dz 



' 2kT r2 



(2) 



where ttlh is the mass of a hydrogen atom. 
The solution is 



where 



P(0) = Pextexp(«;74/i^), 



is the pressure at z = 0, and 



2r ( kT 



1/2 



(3) 



(4) 



(5) 



is the vertical scale height. 

A model for the vertical ionization struc- 
ture is used which is similar to that in Mal- 
oney (1993). The gas is assumed to be highly 
ionized above height Zi where 



rw/2 

L ' 



arecntot{z)dz 



■'i,exi 



(6) 



and 0i,ex is the ionizing flux, which has been 
measured at low redshift by Kulkarni & Fall 
(1993) using the proximity effect. The in- 
ner layer of gas below Zi is shielded and re- 
mains neutral. This sandwich structure oc- 
curs out to a fairly constant column density 
of about 3 X 10^^ cm~^ for each galaxy. For 
r > , where the column density is be- 
low this critical value, ionization equilibrium 
is assumed, where a^ec^tot — C^h for neu- 
tral hydrogen density uhi C is an average 
ionization rate, and arec — 4.07 x 10^^ cm^ 



s-^(T/10, 000K)-°-^5 is the recombination co- 
efficient. 

For lines of sight intersecting the inner 
disk, the observed column density is domi- 
nated by the neutral layer of gas, so that 

where the integral is performed along the line 
of sight (los) for \z\ < Zi. For highly ionized 
gas (r > Tcr) the neutral column density be- 
comes 



P(0) 



C 2kT Jlos 



exp{-2z'^/hl)dz. (8) 



3. Results and Discussion 

The first set of simulations each consisted 
of 20,000 lines of sight through a box con- 
taining 20,000 galaxies. The following stan- 
dard parameters were used for each simula- 
tion, with the exceptions listed in Table 1. 

• The slope of the neutral column density 
distribution was e = 1.5 (1.7 for run 1). 

• The slope of the B scale length distri- 
bution was n — —2 (—1.35 for run 2). 

• The cutoff in the surface brightness dis- 
tribution was at I1Q — 2AB mag arcsec"^ 
(25 for run 3; 27 for run 9) . 

• The ionizing flux 0j_ex = 4.6 x 10^ 
cm~^s~^ (2.8 X 10^ cm~^s~^ for run 4). 

• The ratio of the 21 cm scale length to B 
scale length was /i2i/^B = 3.4 (1.7 for 
run 5). 

• The external pressure was Pext/k — 10 
cm~^ K (1 cm~^ K for run 6). 
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• The switch from exponential to power 
law falloff in total column density oc- 
curred at = ijsyj = 2rcr for run 
7). 

• The luminosity function was changed 
from a Shechter form to a power law 
at Mbsw = —16 (—15 for run 8). 

3.1. Galaxy Counts and the Cosmo- 
logical Density Parameter 

The number density of galaxies for each 
simulation, listed in Table 1, is that which 
is required to explain the observed absorber 
counts of {dN/dz)o = 24.3 ± 6.6 (Bahcall et 
al. 1996) for Nhi > 10^*^-^ cm~^. For galaxies 
with relatively large absorption cross sections 
(run 1) fewer galaxies are required, while the 
opposite effect is seen in runs 4, 5, 6, and 
7. Since the central column density of the 
disk is assumed to vary rather weakly for a 
large change in surface brightness, LSB galax- 
ies have bigger absorption cross sections at a 
given Mb compared to HSB galaxies, so that 
fewer galaxies are required if there are more 
LSB objects (runs 3 and 9). Absorption cross 
sections are generally found to increase with 
Mb, however, so that runs (2 and 8) where the 
shape of the luminosity function is changed 
to favor more luminous objects require fewer 
galaxies. Note that not all dwarf galaxies are 
required to be rich in gas as in these mod- 
els, as the number density of galaxies is rela- 
tively insensitive to the faint end slope of the 
luminosity function. Some dwarfs may have 
formed recently or by varying mechanisms 
(Hunsberger, Charlton, & Zaritsky 1996). 

The simulated luminosity functions (LFs), 
which are normalized to produce absorber 
counts, are shown in Figure |1]. The Shechter 
function found by Marzke et al. (1994) and 



the binned LF from Sprayberry et al. (1997) 
are also shown. The Marzke et al. (1994) 
function generally includes galaxies with /io < 
22 B mag arcsec"^, while the Sprayberry et 
al. (1997) points are restricted to fiQ > 22 B 
mag arcsec"^, so that the sum of the two can 
be considered to be a fairly complete LF. Note 
that only models with relatively small galaxy 
absorption cross sections can provide a rea- 
sonable fit to the observed LF. In other words, 
models where galaxies have larger cross sec- 
tions will overproduce dN/dz if they contain 
as many galaxies as are known to exist ac- 
cording to observed LFs. Therefore, by com- 
paring with observed LFs, it is possible to rule 
out models with overly large galaxy cross sec- 
tions, or put an upper limit on the character- 
istic size of absorbing galaxies. For example, 
a galaxy with = —18.9 would have a ra- 
dius of around 100-300 kpc in runs 4 and 5 
down to a column density of 2 x 10^^ cm~^, 
where some variations occur due to variations 
in central surface brightness. Smaller galaxy 
disks cannot be ruled out by this method, but 
then significant nongalactic absorption must 
also occur. Many parameter assumptions are 
also uncertain. For example, several galaxies 
with /io ~ 27 -B mag arcsec"^ are known to 
exist, so there must be at least some tail at 
the LSB end of the surface brightness distri- 
bution. In reality, some adjustments or vari- 
ations in all of the parameters are expected. 
It can be concluded, however, that it is easily 
possible to obtain a reasonable fit to the ob- 
served luminosity function, and thus explain 
absorbers using galaxies, within the range of 
parameter uncertainties. 

A value for the cosmological density pa- 
rameter VLom, due to the given number den- 
sity of galaxy halos required to match the ob- 
served {dN/dz)Q, is also shown in Table 1. It 
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Fig. 1. — Simulated galaxy luminosity func- 
tions, which arc normalized to produce ob- 
served absorber counts, are plotted along with 
observed LSB (binned points) and HSB (solid 
line) luminosity functions. Note that most 
of the runs shown here require fewer galaxies 
than are known to exist in order to explain 
the observed value of {dN/dz)o. 

is assumed that a fiat rotation curve contin- 
ues to three times the radius of the pressure- 
gravity confinement transition in each galaxy. 
Zaritsky et al. (1997) find that luminous spi- 
ral galaxy halos continue out to at least 200 
kpc, which is consistent with the radius as- 
sumed for a typical L* galaxy in run 4 or 
5. The value of Qdm is generally found to 
be larger for models with a greater number 
density of required galaxies, although varying 
the shape of the faint end of the LF has lit- 
tle efi:ect (runs 2 and 8). Simulations which 
produce enough galaxies compared to the ob- 



served LF tend to have high values of ^dm- 

However, the value of Qdm is most sensitive 
to the shape of the bright end of the LF and 
the properties of the halos, which need to be 
modeled more accurately. 

The value of Qb due to hydrogen in galax- 
ies, which is also shown in Table 1, is more 
closely related to absorber counts. Here Qb 
is shown due to galaxies assuming each one 
cuts off at a neutral column density of 2 x 10^'' 
cm~^, although varying the cutoff column 
density in each run has little effect. Assum- 
ing each galaxy has associated gas out to a 
column density of 10^^ cm^^ increases by 
less than ~ 1%, while assuming that the gas 
cuts off at the pressure-gravity confinement 
transition decreases fi^ by < 10%. Previ- 
ous estimates of Qb due to observed galaxies 
at low redshift have been much smaller than 
those from primordial nucleosynthesis predic- 
tions (Persic & Salucci 1992), but compar- 
ing with the value from Walker et al. (1991), 
^B^ioo — 0.013 ±0.003, most of the runs give 
baryon densities which are somewhat high. 
One exception is run 1, which also requires 
an unreasonably small number of galaxies as 
discussed above. Runs which give particularly 
high Qb include 4, 6 and 7. In run 4, a large 
amount of ionized gas is required, and in run 
7, the gas is too highly concentrated in the 
centers of galaxies. In run 6, the disks must 
be confined by gravity out to larger radii, and 
a more rapid falloff in neutral column den- 
sity occurs in the gravity dominated regime 
(Charlton et al. 1994), so again the gas is too 
centrally concentrated. It can be argued from 
the results of run 7 that galaxies must have 
surrounding gas which falls off slowly with ra- 
dius (such as a power law rather than an expo- 
nential) in order for galaxies to explain most 
of the Lycc absorbers. The only remaining 
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Fig. 2. — Simulated neutral column density 
distributions for low redshift Lya absorption. 
Their shapes are determined primarily by the 
assumed slope, and therefore, and observed 
neutral column density distribution will be 
useful primarily for constraining this slope. 

way to allow for a reasonable fit to the LF is 
to decrease the 21 cm scale lengths as in run 
5, although the value of remains slightly 
high. 

It is possible that the amount of gas in 
galaxy centers has been overestimated by as- 
suming exponential total column density pro- 
files, since the observed central column densi- 
ties tend to be about an order of magnitude 
smaller than the extrapolated fits to exponen- 
tial profiles (Cayatte et al. 1994; de Blok et 
al. 1996). Most of this unobserved gas mass is 
likely to exist in the form of molecular gas or 
stars (McGaugh & de Blok 1997) so that fls is 
not severely overestimated. A somewhat high 



value of flB obtained from the simple mod- 
els used here cannot be taken too seriously 
for a given set of parameters, and small ad- 
justments of multiple parameters can reduce 
the value of Qg. For example, adding more 
LSB galaxies would be reasonable, as some 
tail must exist at the LSB end of the sur- 
face brightness distribution, but the shape of 
this tail is highly uncertain. Also, decreasing 
fsw/fcr might make the models more realistic. 

The number of Lyman limit absorbers is 
shown for each simulation in Tabic 1, again 
assuming that the total number of absorbers 
from Bahcall ct al. (1996) is produced by 
galaxies. The number of Lyman limit ab- 
sorbers has been measured at redshifts > 
0.36, and the lowest redshift values available 
are within the range of {dN/dz)^ ~ 0.2 
to 1.3 (Lanzetta, Wolfe, & Turnshek 1995b; 
Storrie-Lombardi et al. 1994; Stengler-Larrea 
et al. 1995). Storrie-Lombardi et al. (1994) 
find a power law evolution in the number of 
Lyman limit systems which decreases at lower 
redshifts, while the Lanzetta ct al. (1995b) 
and Stengler-Larrea et al. (1995) measure- 
ments are less sharply decreasing with red- 
shift. Assuming that the evolution is flat or 
decreasing to 2; ~ 0, the values of {dN/dz)o^LL 
arc slightly high for most runs. Runs 4, 6, and 
7, which also had unreasonably high fls, have 
excessive numbers of Lyman limit absorbers. 
Thus the same model parameters seem to be 
favored by comparing with Lyman limit ab- 
sorber counts as with as discussed above. 

Since there arc likely to be small variations 
in all of the model parameters, no attempt is 
made to find a single best fitting set of pa- 
rameters. Further discussions concentrate on 
the relatively reasonable run 10, which uses 
the standard parameters except for a surface 
brightness cutoff at 25 B mag arcsec"^, a scale 
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length ratio of /i2i/^B = 1-7, and a change to 
a power law falloff in column density at = 
O.Sr^r- In run 10, 20,000 galaxies and 30,000 
lines of sight are simulated, and a number 
density of 2.0 galaxies Mpc~^ is required to 
explain absorber counts, giving Qdm = 0.37, 
CIb = 0.017, and {dN/dz\LL = 1-56. Note 
that the Lyman limit systems are now pro- 
duced in regions with power law falloffs in 
column density. The number of Lyman limit 
systems would be smaller if the column den- 
sity profiles were modeled more realistically, 
so that the Lyman limit regime is produced 
in a region with column density falloff which 
is faster than the power law, yet slower than 
the exponential. 

3.2. Implications for Lya Absorbers 

Simulated neutral column density distribu- 
tions are shown in Figure ^ (for larger simula- 
tions consisting of 30,000 galaxies and 30,000 
lines of sight in runs and 5). Little variation 
is seen between different simulations since the 
normalization is fixed and the slope is not ex- 
pected to vary much from —1.5. A measure- 
ment of the observed slope at low redshift will 
be interesting, as the value of Qb is sensitive 
to small changes in e, although constraints on 
the other model parameters are unlikely to 
be possible using an observed neutral column 
density distribution. A gap is seen in each 
simulation at ~ 10^^ cm~^ due to the sud- 
den increase in ionization, although it may be 
smoothed out in reality due to dumpiness in 
the gas or variations in the amount of Ionizing 
radiation. 

Absorption cross sections are found to in- 
crease with galaxy luminosity, as shown in 
Figure |^. Thus most absorbers will arise 
from luminous galaxies, although less lumi- 
nous objects may produce some absorption 



3.0 



2.5 - 



2.0 - 

'a 

Q- 

1.5 - 



1.0 - 



-21.0 -20.0 -19.0 -18.0 -17.0 -16.0 -15.0 -14.0 -13.0 -12.0 
Mb 

Fig. 3.— The absorbing (> 2 x 10^^ cm'^) 
radii of galaxies (simulated in run 10) increase 
with galaxy luminosity, although variation oc- 
curs at a given Mb which depends upon sur- 
face brightness. 

since they are more numerous. It was also 
shown by Bowen et al. (1996) that absorp- 
tion cross sections must depend upon galaxy 
luminosity. They argued that it is unlikely 
that low luminosity galaxies produce absorp- 
tion due to this dependence, but a modest 
fraction of absorption at low impact parame- 
ters is produced by galaxies which are low in 
luminosity, as shown in Figures ^ and ^. Vari- 
ations in surface brightness have less effect 
upon absorption cross section, but at a given 
Mb a lower surface brightness galaxy will be 
larger since the column density is assumed to 
vary weakly with surface brightness. Thus the 
upper points in Figure ^ show a sharp edge 
due to the unreahstic sharp cutoff assumed in 
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Fig. 4. — The fraction of absorbing (> 10^^ 
cm~^) galaxies at impact parameters < 300 
kpc is shown function of Mb for run 10. 
LSB galaxies, with ^iq > 22 B mag arcsec"^, 
may cause most of the Lya absorption. Low 
luminosity dwarfs make a modest contribu- 
tion to absorbing galaxies. 



the surface brightness distribution, while the 
lower (HSB) points fall off more smoothly. 

A surprisingly large fraction of absorbers 
(> 10^^ cm~^) is shown to arise in LSB galax- 
ies, as is seen in Figure ^. While this result 
is quite sensitive to the shape of the surface 
brightness distribution and to the relationship 
between column density and surface bright- 
ness, luminous HSB galaxies may make a rel- 
atively unimportant contribution to Lya ab- 
sorption. While LSB galaxies are thought to 
have relatively low central column densities, a 
strong correlation is not seen between column 
densities and surface brightnesses in de Blok 



et al. (1996). Making this relationship weaker 
results in a larger contribution to absorption 
from LSB galaxies. In Figure ^, the absorb- 
ing (> 10^^ cm~^) galaxy impact parameters 
are shown versus Mb- More luminous galax- 
ies can absorb at larger impact parameters, 
as is consistent with Figure |^. It can also be 
seen that low luminosity galaxies make an im- 
portant contribution to absorption if they are 
close to lines of sight and that a large frac- 
tion of absorbers arise in somewhat luminous 
LSB galaxies which allow for absorption at 
the largest impact parameters. Radio surveys 
for absorbing galaxies such as van Gorkom et 
al. (1996) find that absorbers include dwarf 
galaxies at relatively small impact parameters 
and luminous galaxies which are farther from 
lines of sight. 

A plot of the galaxy impact parameters 
versus neutral column density is shown in Fig- 
ure ||. While each galaxy is assumed to have 
a systematic falloff in column density with ra- 
dius, a large amount of scatter is seen in the 
plot mostly due to variations in absorption 
cross sections of galaxies. Variations in disk 
inclinations contribute less to the scatter, so 
that any scenario in which absorption cross 
sections vary would produce a similar plot. 
However, when realistic selection effects are 
introduced, limiting the galaxy sample to a 
more narrow range of properties, less scatter 
is seen. A strong statistical anticorrelation 
is found for all (963) points simulated with 
Nhi > 10^^'^ using a Kendall rank correla- 
tion test. However, when random samples 
of 30 points are chosen, a > 2cr anticorre- 
lation is found in only 122 out of 200 sam- 
ples. Furthermore, those samples which are 
strongly anticorrelated typically contain one 
high column density (> 10^° cm~^) system. 
Lanzetta et al. (1995a) observed an anticor- 
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Fig. 5. — Impact parameters between lines 
of sight and absorbing (> 10^^ cm~^) galax- 
ies in run 10 versus Mb- Luminous galaxies 
can absorb out to larger impact parameters as 
is consistent with Figure 3. Note that most 
objects at low impact parameters are low in 
luminosity. 



relation between equivalent width and corre- 
sponding galaxy impact parameter. Bowen 
et al. (1996) saw no obvious ant icor relation, 
which was used to argue that absorbers are 
more weakly associated with galaxies. Le 
Brun et al. (1996) observed luminous galax- 
ies close to lines of sight, to a lower equivalent 
width limit. They found a weak anticorre- 
lation between equivalent width and impact 
parameter. All of these results appear con- 
sistent with the scenario described here, and 
many more observations will be required to 
distinguish between absorber models. 

The fraction of galaxies which cause ab- 



16.0 



15.0 



14.0 



0.5 




1.0 1.5 2.0 

log(lmpact parameter/kpc) 



Fig. 6. — Impact parameters between absorb- 
ing galaxies and lines of sight in run 10 are 
plotted versus neutral column density for all 
galaxies (circles) and for HSB {^q <22 B mag 
arcsec"^) galaxies (diamonds). 



sorption is shown versus galaxy impact pa- 
rameter in Figures 0a {L > L*, Nhi > W^'^'^ 
cm-2) and 0b (L > O.IL*, Nhi > 10^^ 
cm~^). Bowen et al. (1996) find that their 
observed galaxies typically cause absorption 
out to 300 kpc, while Lanzetta et al. (1995a) 
find a smaller typical absorption cross section 
of about 160 kpc. This value is quite consis- 



tent with cross sections found for HSB, 
galaxies. Lanzetta et al. (1995a) find that 
5/5 galaxies with impact parameters < 70 
kpc cause absorption, where these galaxies 
are slightly less luminous than the ones plot- 
ted in Figure |^. Le Brun et al. (1996) find 
that all of 4 galaxies cause absorption within 
~ 90 kpc where their selection effects are sim- 
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Fig. 7. — The covering factors, or fractions 
of galaxies causing absorption in a given im- 
pact parameter range are shown for run 10. 
The fraction of absorbing LSB (/xq > 22 B 
mag arcsec"^, dashed line) and HSB (/io < 
22 B mag arcsec"^, dotted line) galaxies in a 
given range of impact parameter which cause 
absorption are also shown separately. In Fig. 
7a (7b) the plotted galaxies have luminosities 
> L* (O.IL*) and cause absorption > lO^"^-^ 
cm"^ (10^^ cm~^). LSB galaxies are more 
likely to cause absorption as they have larger 
cross sections at a given luminosity. 

ilar to those in Figure Bowen et al. (1996) 
find that 40 percent of galaxies cause absorp- 
tion at impact parameters between 100 and 
300 kpc, although their sample is not quite 
complete to O.IL* as in Figure ^b. While 
selection effects which depend upon redshift 
may explain some of the discrepancies, it is 
also likely that dwarf and LSB absorbers are 
being misidentified. In the simulation results 
here, the correct absorbing galaxy is known, 
while observers are likely to match an absorp- 
tion line with the most easily visible galaxy. 

The clustering properties of the galaxies 
would need to be included in the models in 
order to test whether galaxy misidentifica- 
tions could help to explain the results from 
Bowen et al. (1996). LSB galaxies may be 



weakly clustered around galaxies which are 
more easily seen, so that the absorption cross 
sections of some luminous galaxies may be 
overestimated. Many galaxies which are more 
difficult to detect may have formed within 
the same extended sheets which contain lu- 
minous, HSB galaxies as seen in cosmolog- 
ical simulations. It is also possible that a 
more clumpy distribution of gas around galax- 
ies, such as that produced by tidal interac- 
tions (Morris & van den Bergh 1994), would 
be more realistic, although the covering fac- 
tor of gas around the galaxies would need 
to remain high out to larger distances in or- 
der to be consistent with observations. Both 
Morris et al. (1993) and Le Brun et al. 
(1996) found that all luminous galaxies within 
0.5 to 1 Mpc cause absorption with a lower 
equivalent width limit. It would be diffi- 
cult to allow for individual galaxies having 
such large cross sections without overproduc- 
ing absorber counts. 

4. Conclusions 

Lya absorber counts at low redshift can 
easily be explained by lines of sight through 
moderately extended galaxy disks when low 
surface brightness galaxies are included. A 
scenario in which absorbers are explained by 
galaxies can easily be made consistent with 
observed luminosity functions, low redshift 
Lyman limit absorber counts, and predictions 
of Qb from standard nucleosynthesis. In suc- 
cessful models the absorption cross sections 
must increase with galaxy luminosity, where 
an L* galaxy typically has an absorbing ra- 
dius of ~ 200 kpc and a low luminosity dwarf 
has a radius ~ 10 kpc. Larger absorption 
cross sections are ruled out by observed lumi- 
nosity functions, assuming that most galaxies 
are potential absorbers. The baryon density. 
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as predicted by primordial nucleosynthesis, is 
overproduced if the gas is too highly ionized 
or too concentrated in the centers of galax- 
ies, while low surface brightness galaxies are 
less centrally concentrated, so that they can 
have larger absorption cross sections without 
containing excessive baryons. 

Low surface brightness galaxies may make 
a large and important contribution to Lja 
absorption. In fact, luminous, high surface 
brightness galaxies may make a relatively unim- 
portant contribution. Luminous galaxies, es- 
pecially those which are low in surface bright- 
ness, are responsible for most Lya absorp- 
tion lines, where the lower column densities 
are often due to luminous galaxies at large 
impact parameters. Low luminosity dwarfs 
are much more numerous, so that their con- 
tribution to absorption becomes more impor- 
tant when they are at smaller impact param- 
eters from quasar lines of sight. Therefore, 
it may be generally difficult to observation- 
ally match a particular absorption line to the 
correct galaxy. 

While a fairly small number of Lya ab- 
sorption lines have been matched to specific 
galaxies to date, these observations gener- 
ally appear consistent with the scenario simu- 
lated here, given that it is impossible to know 
that any particular absorption line has been 
matched to the correct galaxy. A substantial 
fraction of Lya absorption is likely to arise 
from lines of sight through low surface bright- 
ness galaxies which are clustered around more 
easily visible galaxies. These LSB galaxies are 
likely to be misidentified, and improved simu- 
lations which consider the clustering proper- 
ties of these galaxies will be needed to quan- 
tify the effects of possible misidentifications. 
A much larger number of observations will be 
required to test possible models (Sarajedini, 



Green, & Jannuzi 1996), and greater sensi- 
tivity to low surface brightness galaxies near 
quasar lines of sight will be required, even in 
places where luminous, high surface bright- 
ness galaxies have already been detected. Re- 
finements to the models here which allow 
for reasonable clustering properties will allow 
for further constraints on the arrangement of 
neutral and ionized gas relative to galaxies in 
the nearby Universe. 

I am grateful to J. Charlton, G. Bothun, 
C. Churchill, R. Ciardullo, S. McGaugh, J. 
Miralda-Escude, and D. Schneider for valu- 
able discussions and to D. Sprayberry for sup- 
plying data points. This work was supported 
by a fellowship from the NASA Graduate Stu- 
dent Researchers Program. 
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Table 1 
Simulation Results 



run # 


parameters varied 


galaxies (Mpc) ^ 






{dN/dz)o,LL 





all standard 


0.62 


0.21 


0.021 


1.59 


1 


e = 1.7 


0.16 


0.06 


0.006 


0.415 


2 


n = -1.35 


0.28 


0.19 


0.020 


1.39 


3 


IIq = 25 cutoff 


0.52 


0.16 


0.018 


1.31 


4 


(/.i,ex = 2.8 X 10^ 


1.4 


0.50 


0.049 


2.25 


5 


h2i/hB = 1.7 


2.8 


0.53 


0.023 


1.41 


6 


Pext/k = 1 


1.9 


1.00 


0.063 


4.77 


7 




2.3 


0.44 


0.077 


4.04 


8 


Mbsw = -15 


0.44 


0.21 


0.021 


1.34 


9 


110 — 27 cutoff 


0.39 


0.11 


0.014 


1.22 



Note. — Columns are: (3) the number density of galaxies required to explain absorber counts, 
(4) the cosmological density parameter due to dark matter halos (out to three times the gravity 
confinement radius in each galaxy), (5) the cosmological density parameter due to hydrogen (out 
to a column density of 10^^'^ cm~^), and (6) the number of Lyman limit systems per unit redshift. 



